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1 Introduction 



In the past few years supersymmetric quantum mechanics based on shape 
invariance and intertwining concepts has manifested a significant progress 
j|l|]-[10|. Its technique started to influence not only the traditional branches 



of physics such as atomic, nuclear and high energy physics, which originally 



stimulated its emergence [11 , but also the classical areas of mathematical 



physics and the theory of differential equations. Recently, in fl2|-|[13|| the 



detailed investigation of factorization technique has been performed for one 
specific form of second order differential equations (SODE) with polynomial 
coefficients, admitting polynomial solutions based on known Rodrigue's for- 
mula [14]]. We choose similar initial arguments to construct explicitly a wide 



class of QM potentials for ID Schrodinger equations admitting (after sepa- 
ration of asymptotic behaviour of the wave function) polynomial solutions. 

Our consideration and analysis in sections 2 and 3 shares a common subject 



with the known Natanzon papers |T5|,|T6[ but differs from them by being 



more general because of the fact that we work with Papperitz rather than 
hypergeometric equation with arbitrary positions of the singular points (in 
complex domain). Moreover, in addition to transformations described in 
[P^[],||16| where the spectral parameters are preserved, we include into consid- 
eration the cases when we change the roles of spectral parameters in original 
polynomial system and in the Schrodinger equation (generalizing in this way 
the known consideration of the Coulomb case, see the discussion in the text). 



One more thing we would like to mention is that Turbiner's approach |T7| 
to generalized Bocher problem is very near to ours for exactly solvable case 
(though in |L7| the possibility of the transformation to appropriate Sturm- 
Liouville problem is merely mentioned rather than investigated in full de- 



tails). But for quasi-solvable cases our approach differs from those in JT7 
because we do not use the factorization ideas and investigate symmetry prop- 
erties but concentrate on the connection of the appropriate polynomials with 
the corresponding family of Schrodinger equations. It allows us to explicitly 
define the additional relations for the operator to be zero-grading (in terms 
of JH]) and gives the way of regular construction of quasi-solvable poten- 
tials (all inside some definite family) with an arbitrary chosen number of 
algebraically constructed eigenstates. 

The last remark we should make is that we do not intend to perform the 
comparison of the proposed approach with all other known methods of the 
construction of exactly solvable quantum potentials as this should be defi- 
nitely the topic of a review paper rather than original research paper. But, 
in our opinion, the proposed approach will definitely influence the reviewing 
of the results obtained in other ways (especially for quasi-solvable cases) and 
further establishment of its inherent relations with other methods. 
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The paper is organized as follows. In the Section 2 we start with the special 
case leading to polynomial solutions, namely the polynomial family intro- 



duced in |]T2|-JT3|. Section 3 is devoted to the explicit construction of the 
Schrodinger equation corresponding to that polynomial family, and to the 
presentation of some classification scheme as well as the discussion of its 
relation to the known solvable cases. In Section 3 we represent the gener- 
alization of the proposed scheme in some directions, namely the application 
to the partial differential equations and to higher order polynomial coeffi- 
cients, and we demonstrate one non-trivial irreducible example and make 
some concluding remarks on the applicability of the considered approach. 



One Construction of Polynomial Solutions 
for SODE 



Practically all solvable ID quantum problems correspond to those which 
can be transformed to the equation of hypergeometric type, and that gives 
eigenvalues of bound states by the requirement of finite hypergeometric series, 



thus being a polynomial of a given order |14|, [18] . We start the consideration 
from the case leading precisely to that known situation, though from a little 
bit another point of view. Recently |l2|-|jl3l it was shown that the eigenvalue 
problem for the operator of the form 

^wMi( A(l) "' ll) S W 

leads to polynomial solutions with special requirements for the functions 
A, W. Namely, if we choose A(x) as a polynomial of at most a second order, 
let us define, 

A(x) = a + a\X + a 2 x 2 , (2) 

and W(x) as a non- negative function such that pp^y^: (A(x)W(x)) is at most 
a first order polynomial 

B(x) = b + bix, (3) 

then we can construct an orthogonal polynomial family being a solution for 
the eigenvalue problem, namely for the operator £, 

A(x)W(x)— + 7n$n(x) = 0. (4) 



W{x) dx \ dx 
The polynomials given by the classical Rodrigue's formula |TJ 
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are orthogonal with respect to the weight function W(x) on the interval (a, b), 
chosen such that the following conditions hold 

A(a)W(a) = A(b)W(b) = 0. (6) 



This can be shown by classical consideration as exposed e.g. in |]14| . In this 
way we can choose the interval as inside the roots of the polynomial A(x), 
if the latter one possesses the real roots, or the roots of W(x), including 
infinity points, in the case of the functions tending to zero at infinity, or 
some combination of both alternatives. The eigenvalues 7„ turn out to be 
given as fI~2]l- [|T3[ 



f (A(x)W(x)y \ n(n-l) , 
In = -n y w{x) j - 2 ^ W 

which in case (0) and (j^) is equal to expression (11). 

At that point, having enough information, we can write and solve the equa- 
tion for W(x), namely 

— (A(x)W(x)) - B(x)W(x) = (8) 

(JjJb 

As the equation is linear ODE of the first order its solution has the form, 
explicitly 

W(x) = exp / = Pl ;^ ; , (9) 

A[x) [J A[x) J ao + a\x + a^x 1 

and, of course, the integral we wrote is easy to calculate 



b + b\x b\ log(a + a\x + a^x 



2\ 



dx = 1 taV u ' 1 ' z ' - (10) 



ao + a\x + d2X 2 2a 
(— 2a2&o + ctibi) arctan( 
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ai+2(i2X 
\/ _a l 2 + 4a a 2 ' 



a2\/ — &i 2 + 4aoct2 

but for our purposes it is convenient to use it in the form we represent W(x) 
in equation @. 

It is evident from the direct substitution that in the case (0) and (|3]) the 
eigenvalues 7 n are given explicitly by 

7n = -n(bi + a 2 {n- 1)), (11) 

while the equation for the polynomials becomes 

(a + a x x + a 2 x 2 )$"(x) + (b + bix)& n (x) (12) 
- n(&! + a 2 (n - l))$ n (x) = 
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The equation (|T2| ) is SODE with 3 singular points at the roots of A{x) (we 
denote them as Xi,X2) and infinity. It is easy to check that if the roots of 
A (x) are different the equation is of hypergeometric type (generally speaking 



Papperitz equation ||14|), whilst it is confluent hypergeometric when the roots 
are coincident. 

The question naturally appearing is what sort of quantum mechanical prob- 
lems could be associated with the polynomial family we described. To an- 
swer it, there is a natural way, namely we can try to implement adjusted pair 
of coordinate transformation and similarity transformation for the equation 
( |12|) in such a way, as to obtain the constant coefficient at second deriva- 
tive and zero coefficient at first derivative. The resulting equation will be of 
Schrodinger type. Let us perform this program. 



3 Transformation to the Schrodinger equa- 
tion 

Let us forget for a while that A(x) and B(x) are polynomials. We have the 
equation 

A 2 (x)y"(x) + A^y'ix) + ey{x) = (13) 

with arbitrary coefficient functions A\(x), A 2 (x). First, we make the variable 
change 

x = F(u), 
d 1 d 



dx F'(u) du 
d 2 1 d 2 F"(u) d 



dx 2 F' 2 (u)du 2 F' 3 (u)du 



(14) 



and choose the transformation in a form allowing to introduce some yet 
undefined but prescribed function of new coordinate oj(u) (which we could 
define later for the sake of the most convenient choice) 

u 2 (u)[F\u)} 2 = A 2 (x) (15) 

to get 

9/ x »/ x »/ x , s (2A 1 (F(u)) - A'JFiu)) \ 
u 2 (u)y"(u) + y'{u)uj{u) U ) U= -- ^ + uj'(u) (16) 

+iy(u) = o 

where prime means the differentiation with respect to the function's argu- 
ment. Now we implement the similarity transformation Y(u) = exp(x{u))y(u) 
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and choose the function xi u ) i n such a way as to kill the term with the first 
derivative, so that we must have 



A 2 \F{u)) Ax(F(u)) 



Mu) \2yjA 2 (F(u)) ^A 2 (F{u)) 
Then, the equation is transformed into 



uj\u) 



u 2 (u)Y"(u) + Y(u) 
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A 1 (F(u)f , A 1 (F(u))A 2 '(F(u)) 



4A 2 (F(u)) 



+ 



2A 2 (F(u)) 



(17) 



(18) 



3A 2 '(F(u)f uo\uf A 2 "(F(u)) A x '{F{u)) uj(u)u"(uY 



+ 



16 A 2 (F(u)) 4 4 2 2 

which can be considered as Schrodinger type equation if we manage to iden- 
tify and to separate some free constant parameter in it at Y(u) (playing the 
role of "energy"), after the division of both terms by u 2 (u). As we see in the 
trivial choice of u(u) = 1 we simply get the Schrodinger equation with the 
energy 7 and the potential V(u) given by 



V(u) 



4 Ai(F{u)) 2 - 8 A\(F{u)) A' 2 (F(u)) - 3 A' 2 {F{u)f 
16 A 2 (F(u)) 

A[(F(u)) A»(F(u)) 



(19) 



2 4 

One known case when we have to choose u(x) not equal to unity is the 
Coulomb potential as we will see below. Now we can use the fact that both 
A 2 (x),Ai(x) are polynomials, choosing them in accordance with (0,0) as 
A 2 (x) = a 2 x 2 + a\x + a and Ax(x) = b\x + b to obtain the equation 

~a 2 - 61 + 2 7 lo'(u) 2 uo"(u) 



Y"(u) + Y(u) 



2uo 2 {u) 

3( ai + 2a 2 F(u)) 



4a; 2 (m) 2w(u) 
2 



+ 



16lu 2 (u) (a + ai F(u) + a 2 F(uf 
a 1 + 2a 2 F(u)) (bo + hFiu)) 



(20) 



2uj 2 {u) (a + at F{u) + a 2 F(u) 

(bp + hFju)) 2 
Au 2 (u) (a + ai F(u) + a 2 F(u) 



Now, we start to analyze first systematically the simplest case putting func- 
tion u(u) to be unity. Then, for the potential of the Schrodinger equation 
we have 

\2 



3( ai + 2a 2 F(u)Y 



V(u) = 1 1 — 

2 2 16 (a + aiF(w) + a 2 F{u) 



+ 



(ai + 2a 2 F{u)) % + 



+ 



(bo + hFiu))' 



2 (a + Oi F(u) + a 2 F(m) 2 ) 4 (a + a x F{u) + a 2 F(u 



(21) 
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We can explicitly find the dependence x = F{u) by solving the equation fll5|) . 
Then, taking the inverse function, we get for x 



ai sinh( v^p" ) + v /ao02 sinh(^/a^u) ^ q -p ^ Q 



a-1 

2a~2 ' ° 2 + U ' ^ - U (22) 

v^M+^f^, a 2 = 0, ax^O 

|= it, a2 = 0, ai = 0, ao 7^ 



where P = a\ — 4aoa 2 is the discriminant for A 2 {x). 

The explicit expression for the quantum potential can be obtained after the 
substitution of (p2|) into expression (^) but for the general case the resulting 



formula becomes too complicated. Specifying the values of the parameters 
a,, 6j it is possible to show that almost all known solvable cases in quan- 
tum mechanics except the Coulomb potential, are inside the potential family 
we constructed. The latter one will be analyzed later when we try to con- 
struct new potentials by choosing a nontrivial u(u) function. So let us first 
investigate the case uj{u) = 1. 

We shall classify the different cases by the roots of the polynomial A 2 (x). We 
have two topologically different cases when A 2 (x) is not degenerate, namely 
when A 2 (x) has two different roots and the discriminant P / 0, and the 
case of one root with degeneracy 2, T> = 0. As one can easily see in the 
space of the parameters ai the first case fills the region inside and outside the 
conical surface for which the equation is T> = 0. So we call the appropriate 
cases regular and irregular (for D = 0), respectively. We shall refer to them 
sometimes as Jacobi and Morse cases, based on the name of the appropriate 
polynomials (for the first one) and solutions (second one). 

The additional cases appear in (^) as a result of the degeneracy of the 
polynomial A 2 (x), it could be of the first order (we call this case Laguerre's 
case; a 2 = 0) and of the zero order (further referred to as Hermite case; 
a 2 = a 1 = 0). 

The first remark we would like to make is that, our polynomial family has five 
parameters whereas spectra depend on top power coefficients of A 2 (x), A\{x) 
only, namely 02,61 (see equation (|lT|)), and so we have the evident freedom 
of choosing some parameters without loosing characteristic features of the 
problem. Obviously, we can change parameters bo, a\ simply by the trivial 
change of the origin in x variable and the scale on it. Then, if we choose 
definite values of a 2 , bi (one of them could be considered as a scale for energy 
and could be chosen e.g. as unity) we obtain a two parametric family of 
polynomials and one parametric family with full isospectrality property. The 
variation of the parameter a$ then will lead to different non-trivial cases we 
mentioned. 
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So, we start with the regular case D / 0. Then the change of variables 
x — > u is given by the top line in (|22"|), and the orthogonal polynomials have 
the weight function given by 

W(x) = P{ a2V 11 )] (23) 

(a + a\ x + a 2 x 2 ) 2a 2 

As this family of polynomials has no commonly used name we will refer to 
it as to generalized Jacobi polynomials, ordinary Jacobi case corresponds to 
ao = 1, ai = 0, a 2 = — l,&o + bi = 2p, b — b\ = 2q corresponding to symmet- 
rically chosen real roots at ±1 and the interval of orthogonality [—1, 1]. 

The quantum potential has the following general form [z = JIiqu) 



■ i A + B sinh z + C cosh z + D sinh 2z + E cosh 2z , . 

V{u) = -2 (24) 

(2^/a a 2 cosh z + a-y sinh z) ) 

with the coefficients A, B, C, D, E expressed through the original ones as 

A = a 2 2 (5 ai 2 - 20 a a 2 + 8 b 2 ) + (3 a x 2 - 4 a a 2 ) Oi 2 

+ 2 a 2 (—3 ai 2 + 12 a a 2 — 4 ai bo) b\, 
B = 8 ^ao a 2 (61 - 2 a 2 ) (2 a 2 6 - ^i h) , 

C = -4 ai (2 a 2 - fo x ) (2 a 2 6 - a a h) , (25) 
£) = 4ai y/a a 2 (a 2 - &1) 2 , 
E = (ai 2 + 4 a a 2 ) (a 2 - i>i) 2 . 



It is straightforward to see that the potential family ( f24|) includes Poschl- 
Teller potentials (both ordinary and modified), Scarf-like potentials, Rosen- 
Morse and Manning-Rosen potentials [[Lj|] at appropriate choice of the pa- 
rameters. 

In singular case D = it is more convenient to introduce other parameters 
rather than a^, namely ao = a 2 , a 2 = /3 2 , a% = 2a/?, automatically satisfying 
the degeneracy condition, and then the weight function becomes 

W {x) = {a + M^e W {^^}, (26) 

and the potential reads through newly introduced coefficient A, B, C 

V(u) = A + Be~ Pu + Ce~ Wu , 



A 



4/3 



2 



(b ia -b (3) (h-2(3 2 ) 

B = 2^ (27) 

_ (ha - b (3) 2 
4a 2 f3 2 ' 
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This evidently corresponds to the Morse class of potentials fll8| 



The case when A 2 (x) becomes the first order polynomial (a 2 = 0), gives for 
the weight the following formula 

W(x) =e^T (ao + dixf ~~^T ? (28) 

and for the potential 

V{u) =A + Bu~ 2 + Cu 2 , (29) 

with the new parameters A,B,C (not to be confused with those obtained 
above) expressed through the old ones as 

A _ b 1 (a 1 b -a b 1 ) 

B = , (dUJ 

U _ 16' 



The resulting potentials, as we see, are the combination of harmonic oscil- 
lator (HO) potential plus centrifugal potential B/u 2 . And for the sake of 
completeness it is worthwhile to mention, that the case when A 2 (x) is a con- 
stant (a 2 = a,\ = 0) corresponds to ordinary HO case with the oscillator 
position shifted by — &o/&i- 

Before going further, let us construct the explicit representation for the wave 
functions of the appropriate Schrodinger equation and let us discuss the 
bound states within this approach. As we made two subsequent transforma- 
tions to obtain the Schrodinger equation, the solution in terms of polynomials 
has the form 



Y(u)=exp{ X (u)}P n (F(u)) (31) 

where xi u ) is given by ( |PTD and F(u) is given by (p~5|) . The energy corre- 
sponding to this eigenfunction turns out to be the sum of 7„ ( see equation 
(|TTD ) and some constant factor depending upon the parameters Oj, 6, and 
leading to the shift of the energy's origin. As the family of orthogonal poly- 
nomials has infinity and countable number of members, though the class 
of the potentials includes not only those which grow indefinitely at infinity 
(supporting bound states only), but also such with a finite number of levels 
(and finite ionization energy), we have to understand what is the condition 
for a bound state in the system. Indeed, this is very simple in the discussed 
case, the function W(u) gives the asymptotic behaviour of the ground state 
wave function and as the point transformation x = F(u) could be non-trivial, 
the resulting high order polynomials P n (F(u)) can have growing behaviour 
at infinity which might more than compensate that of W(u) and thus makes 
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Y(u)'s norm infinite. Therefore, the condition for the bound states is simply 
the condition of a finite norm of Y(u). Of course, an interesting question ap- 
pears whether the non-normalizable solutions of polynomial type correspond 
to some physically significant features of the system, e.g. to quasi-bound 
states (long-lived localized states) embedded in the continuum, but we will 
not discuss it here. 

Now, we can consider other possible choices of the function oj(u). This is 
stimulated by the known sequence of transformations for the Coulomb prob- 
lem [|TS|], where the first step is the change of scale in a way to put "energy" 
parameter into potential function with subsequent transformation of the orig- 
inal Schrodinger equation into hypergeometric equation. 

As we can see from the equation ( p0|) the problem is pure algebraic and 
there are several ways to try to obtain the free parameter which could be 
interpreted as energy. The first one is to choose u{u) to satisfy the equations 
either u'{u)/u{u) = k or uj"(u)/uj(u) = k, that will produce constant factor 
due to fractions including the appropriate ratios in equation (^t]). The second 
one takes place when uj 2 {u) = A 2 (F{u))~ k , k > 0, that could lead to free 
coefficient in potential due to cancellation of some denominators in (]2"0"D. We 
shall not pursue these cases any further, but treat the most important case 
instead. Namely, if the order of polynomial A 2 (x) is less than two, then new 
possible cases also appear, as we will see, for k = —1, which turns out to be 
precisely the Coulomb case that we discuss below. The last and more special 
case is realized when A 2 (x) has different real roots and the coefficients bo, b\ 
are chosen in such a way as to construct common divisor (of the first order) for 
both numerators and denominators in potential. In this case choosing u(u) 
in the form of u(u) = (F(u) — x±) we also obtain free coefficient in potential. 
We will not consider all the above mentioned cases in detail here but restrict 
ourselves to one specific choice of A 2 (x), stimulated by the Coulomb problem. 

Let us assume that A 2 (x) = x and we will choose the uj(u) in the form uj 2 {u) = 
F(u)~ k . Then, if we choose k = —1, it is easy to see from the equation (|15|) 
that the point canonical transformation turns out to be identity, and we get 
the standard Coulomb case 



The case k — lis also a special case here, so we have the following expression 
for the variable change x = e~ u and the Schrodinger equation takes the form 




boh 



2u 



+ 1 Y{u) = 0. 



u 



(32) 



Y"(u) + 



(bo ~ I) 2 



(6o6i-2 7 ) 



— u 




-2u 



Y{u) = 



(33) 



4 



which is the case of the quantum Morse potential. 
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For k ^ ±1 we obtain after integration of the equation fll5|) the following 
expression for the point transformation, 



x 



(fc-iXit-d) 



2 



(34) 



The substitution of the last expression leads to a fairly complicated form 
of the Schrodinger equation for Y(u), but as one can show, there are no 
more cases except those we mentioned, where it is possible to obtain a free 
parameter in the role of " energy" . 

The successful implementation of the construction of solvable Schrodinger 
potentials, as we already have seen, was due to the evident existence of poly- 
nomial solution of the equation ([121) . It is possible to find the generalization 
in more complicated cases, which is the topic of the following section. 



4 Some generalizations of the approach 

As we saw, the main feature of the system considered was that the second 
order differential operator £ (see equations ([!]) and (|T2"D) preserves the linear 
subspace M. n of the polynomials of order n for all n. It was due to the 
special adjustment of the orders of polynomial coefficients with the order of 
appropriate derivatives. This idea, of course, can be implemented not only 
for a special case of SODE and polynomial coefficients up to the second order, 
but in much more general case of linear PDEs. Indeed, we can construct the 
following general form of the n-th order linear differential operator £ with 
the same property, so that it preserves the space when acting in the space 
of the polynomials of m variables x = {x±, ...x m }. The general form of the 
appropriate linear PDEs reads 

n 

£Y{x) = "£P j+N {x)d K .Y(x) = 0, (35) 

where N is some non-positive integer number (when iV < we have de- 
generate cases, see analogous discussion on SODE in the previous section). 
We introduce multi-index Kj cis J- th order partial derivative over arbitrary 
combinations of variables in a standard way, by 

We also say that the weight of Kj equals j and write it as jj^Kj = j. It is 
worthwhile to point out that we even do not need to demand the commu- 
tativity of the derivatives, so that the same consideration could be applied 
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for non-commutative case (quantum groups and quantum algebras see, e.g. 
[pCp ). Moreover, we can consider the operators which do not preserve such 
finite dimensional spaces, but map one into another (with higher dimension). 
The latter is just the permission for N to be positive. 

In the case of standard partial derivatives, when we are looking for the so- 
lution in terms of k-th order polynomials in the ring F[x\,X2, ■-,x m ], every 
term in the sum in (|55| ) maps the argument into the space spanned by the 
monomials x KN+k = i'/x^ 2 ...i m tm , {h + • • • i m = N + k}. The latter space 
is a finite dimensional vector space and a direct sum of the spaces of sym- 
metric homogeneous polynomials corresponding to different permutations of 
indices for the monomials written above. We denote the space spanned by 
the definite monomials of order k as ^ k '%i 1 ..,i m ]) ■ Then we can write down the 
expression for the dimension for the image space for operator action 

k+N m 

dimM%=Y. E dvm^r [h .„ im] . (37) 

ilH \-im=m 

Now, the construction of the k-th order polynomial solution 

Y k {x)=j^C Kj x^ (38) 

3=0 

leads simply to the linear algebraic problem for non-trivial solution for the 
coefficients C Kj . 

At this point two different cases are possible. The first one is realized if 
N < 0, that is the maximal order of polynomial coefficient is less or equal to 
the PDE's order. In this case we can always satisfy the system of equations 
because the number of linear homogeneous equations for C K] is precisely 
equal to dim -M°. Then, the non-triviality condition is the condition of zero 
determinant for corresponding matrix obtained from equating all coefficients 
at monomials of type x Kl to zero, and this gives us spectral parameter for 
the polynomial family, namely the quantization condition imposed on the 
coefficient Pq(x). Then, considering the problem over the field of complex 
numbers, we can always construct a polynomial family in this case for some 
quantized value of the coefficient Po{x). In the contrary, when we have the 
condition N > 0, we are still obliged to fulfill dim M. k conditions but only for 
dim M° k coefficients C Kj . The resulting system becomes overcomplete which 
simply means that we can construct some separate polynomial solutions of 
the equation ( |35] ) for only a few levels, maybe even one, that is for definite 
choice of n and, additionally for special values of some of the coefficients in 
the coefficient polynomials. 

It is very interesting to mention here that in the 1-D case considered in the 
previous section for SODE, the appropriate matrix turns out to be upper- 
three-diagonal, with additional relationship between elements, so that its 
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determinant for j-th order system has the following form (preserving notions 
for coefficient of A 2 (x), A%(x)) 

( 7 b 2a ... 

7 + 61 26 ... 

••• (j-l)6 + (j-l)(j-2)a 1 j(j-l)a 

... 7 + (j-l)6 1 + (j-l)(j-2)a 2 jb +j(j-l)ai 

V ... l + -l)a 2 

that leads to one degenerate eigenvalue (rather than j + 1) for a given j-th 
order polynomial. 

What realizations of the scheme described above could be successfully used 
for the construction of the solutions for Schrodinger equation except those 
which we demonstrated in Section 2? There are three evident but not easy 
ways. The first and the easiest one is the consideration of ID problems with 
higher order polynomials, and we intend to demonstrate this in one example, 
putting away its full description as a subject of separate publication. We 
consider the construction of the polynomial solutions starting from the third 
order polynomials to give the representation of the problems emerging there. 
Let us have the equation of the form 

x 3 y"(x) + a(x 2 - l)y(x) + {J3x + i)y(x) = 0. (39) 

Then, in the same manner as we did in Section 2 we map the equation to the 
Schrodinger equation via two subsequent point transformations x = 4/u 2 , 
and the gauge transformation 

X (u) = Ci + — + — ^— log u. (40) 
The resulting equation has the form 

Y"(u) + (7 - V (it)) Y(u) = 0, (41) 

rr , x (3a-a 2 )u 2 u 6 a 2 4 a 2 - 8 a - 16 (3 + 3 . . 

V(u) = 1 1 = . (42) 

v ; 8 256 4m 2 v ; 



The solution will be given by the formula (|3lD, but now we have to construct 
the polynomials in a non-trivial way, because the Rodrigue's formula is no 
more applicable. So we start to look directly for the polynomial solution of 
the equation ([39]) . Let us define n-th order polynomial as p n (x) = Y^I=q c % x% ■ 
When we substitute this anzatz into equation (j39"D and equate to zero the 
coefficients at all orders of independent variable x, we obtain n+1 equations 
for n coefficients q, i = . . . n — 1, 7 (the last coefficient c n should be chosen 
in order to satisfy the standardization condition). Therefore the system is 
overcomplete and for a non-trivial solution we must specify some additional 
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coefficient in a unique way or to impose one additional condition on some 
parameters of the system. In our case we have the only choice to add f3 to 
the list of coefficients to be found. Let us find then the solution for example 
for n — 1, 2 explicitly starting from the case n — 1. 

Solving the set of equations for c ,(3, 7, two solutions Cq = ±1, (3 = 
—a, 7 ± = i=a can be constructed that giving for Yi(u) 



Y±(u) = e^ (zpl + ^ju^r 2 , (7 = T a). (43) 



We have to impose the finiteness condition for the norm of the solution, so 
that in our case a < —1/4. The last requirement follows from simultaneous 
demand on proper behaviour at infinity that leads to a < and integrability 
at u — 0. But as one can see, in the interval —1/2 < a < — 1/4 the potential 
becomes repulsive as u — > 0. The last property means that we have to impose 
the boundary condition Y(u)\ u=0 = that leads to the restriction a < —1/2 
on the admittable region for the variation of a. 

Then, both solutions will be, as one can see, bound states of the system and 
the solution for larger 7 will have zero not only at u = but also at u — 2, 
which evidently corresponds to the first excited state, whereas the first one 
corresponds to the ground state. 

In a similar way, for n = 2 we obtain the condition f3 = — 2 (1 + a) and the 
following three solutions for c , c±, 7 

4 ] = -^T, c? = 7 (1) = 0, (44) 
a + 1 

ct + 2 a + 2 v 

Then, the appropriate eigenfunctions are given by 

y»(„) = e *„i-(>«- _^ ), (45) 
\m 4 a + 1/ 




.,„• , ,u 4,/2a(2a + 3) 

( u ) = e — «2 " a — + ± 1 



a + 2 (a + 2)m 2 



The similar consideration as we performed for n = 1 shows that admittable 
region for the parameter a now is given by a < —5/2. Then the constructed 
eigenstates represent the ground and the first two excited states for the po- 
tential 

Tr , x a 2 u 6 a(a-3)u 2 4a 2 + 24 a + 35 

V{U) = "256" - 8 + 4^ ' < 46) 

The remarkable feature of the example we considered was that we construct 
some eigenstates corresponding to a given potential using polynomial anzatz 
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of a given order. The solution for those cases corresponding to the ground 
state and to the excited states was given by the number of the roots of the 

polynomials of a given order. This is in contradistinction with the standard 
solvable situation where polynomial's order is equal to the quantum number 
because for the classical polynomial family, the n-th order polynomial has 
precisely n roots on an interval where the family is defined, unlike in the 
general case, where the real polynomial can have complex (nonreal) roots. 

^From the last we can make some useful conclusions. As we see, if the 
order of polynomial coefficient functions is different from the order of the 
equation, the polynomial solutions in general do not exist except for the 
special values of the parameters. In the example that has been considered, 
we have to put one additional condition on the parameter to construct non- 
trivial solution. Nevertheless, after appropriate restriction of the region for 
some parameters (a in the considered example) and by fixing the value for 
some other (expressing j3 through the a in the discussed case) we were able 
to construct eigenfunctions of some low lying eigenstates 

Now, we return to the discussion of other possible generalizations of the 
proposed approach. 

The second admissible choice is to consider, in the same way as we did, the 
Schrodinger equations solvable in momentum representation. As it is easy 
to understand, physically interesting potentials of polynomial type with the 
order greater than two will correspond to higher order differential equations, 
so that we can e.g. ask which Schrodinger equations can be constructed, 
based on the polynomial solutions for the equation like this A±{x)y"" \x) + 
A 3 (x)y"' + ... = 0? We stop the discussion of this possibility at this point 
leaving it for the future publications too. 

The last and the most complicated case corresponds to the consideration of 
genuine higher dimensional problems and appropriate PDEs. The most diffi- 
cult obstacle to be overcome there is the necessity to perform transformations 
from original equation for polynomials to the equation of the Schrodinger 
type. Although the theory of characteristics is applicable in this case, the 
resulting equation, at first glance, could be hardly interpreted in terms of 
the Schrodinger type. As for the latter one we must demand the existence of 
pure constant term included into coefficient function at zero-derivative (see 
the discussion after equation (pT8[)). Nevertheless, this direction is of great 
importance for deeper insight into integrability and solvability problems in 
quantum mechanics. 
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5 Discussion 



In summary, we demonstrated that one can construct explicit formulae for 
the family of the orthogonal polynomials depending on five parameters, and 
thus we can associate with them the family of isospectral potentials (isospec- 
trality with respect to 3 free parameters) which include almost all known 
quantum mechanically exactly solvable potentials. Some generalization of 
the approach to higher dimensional equations (PDEs) as well as to the higher 
order ODEs has been proposed. 

We may conclude with a little bit speculative thought. If one would put for- 
ward the requirement for the bound states of a quantum system (in analytical 
case) as a demand on polynomial type of the reduction of the wave function 
(that seems to be reasonable and evidently fulfilled for all up-to-date known 
ID solvable cases), then the immediate conclusion follows that the proposed 
approach (with its generalizations discussed) includes all analytically solvable 
cases. 
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